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ABSTRACT 

We present panoramic Spitzer MIPS 24-/im observations covering ^9x9Mpc (25' x 25') fields around two 
massive clusters, CI 0024+16 and MS 0451-03, at z = 0.39 and z = 0.55 respectively, reaching a 5-cr flux limit of 
^ 200/iJy. Our observations cover a very wide range of environments within these clusters, from high-density 
regions around the cores out to the turn-around radius. Cross-correlating the mid-infrared catalogs with deep op- 
tical and near-infrared imaging of these fields, we investigate the optical/near-infrared colors of the mid-infrared 
sources. We find excesses of mid-infrared sources with optical/near-infrared colors expected of cluster members 
in the two clusters and test this selection using spectroscopically confirmed 24/im members. The much more sig- 
nificant excess is associated with CI 0024+16, whereas MS 0451-03 has comparatively few mid-infrared sources. 
The mid-infrared galaxy population in CI 0024+16 appears to be associated with dusty star-forming galaxies (typ- 
ically redder than the general cluster population by up to Ay ^ 1-2 mags) rather than emission from dusty tori 
around active galactic nuclei (AGN) in early-type hosts. We compare the star-formation rates derived from the 
total infrared (8-1000/im) luminosities for the mid-infrared sources in CI 0024+16 with those estimated from a 
published Ha survey, finding rates > 5 x than those found from Ha, indicating significant obscured activity in the 
cluster population. Compared to previous mid-infrared surveys of clusters from z ^ 0-0.5, we find evidence for 
strong evolution of the level of dust-obscured star-formation in dense environments to z = 0.5, analogous to the rise 
in fraction of optically-selected star-forming galaxies seen in clusters and the field out to similar redshifts. How- 
ever, there are clearly significant cluster-to-cluster variations in the populations of mid-infrared sources, probably 
reflecting differences in the intracluster media and recent dynamical evolution of these systems. 

Subject headings: galaxies: clusters - galaxies: starburst - galaxies: infrared - clusters: individual CI 0024+ 16, 
MS 045 1-03 
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1. INTRODUCTION 

The galaxy populations within the virialised regions of rich 
clusters at z are characterised by passive elliptical and 
lenticular (SO) galaxies (Oemler 1974; Dressier 1980). In con- 
trast, 5-Gyrs ago, at z 0.5, the galaxy populations in the most 
massive clusters had larger fractions of star-forming late-type 
spirals, and a corresponding deficit of luminous SO galaxies 
(e.g. Dressier et al. 1997; Small et al. 1997; Couch et al. 1998; 
Fasano et al. 2000; Treu et al. 2003). Taken together, these two 
observations imply that a process (or processes) is transform- 
ing many of the star-forming, late-type spirals in these regions 
into the passive early-type population (specifically SOs) found 
in local clusters (Poggianti et al. 1999; Moran et al. 2006, in 
prep). 

When considering potential pathways to produce this evo- 
lutionary change, we need to bear in mind that the typical lu- 
minosities of the star-forming spirals appear to be too low for 
them to transform into typical SO galaxies found in local clus- 
ters, without the addition of significant numbers of new stars 
(Poggianti et al. 1999; Kodama & Small 2001). This problem 
is exacerbated when we include the fading which is likely to 
take place after the cessation of star formation in these galax- 
ies. This then leads us to concentrate on mechanisms which are 



capable of increasing the luminosity of the galaxies - mergers 
and starbursts. There has been a recent upsurge in interest in 
the potential for so-called "dry" mergers (mergers between dis- 
sipationless stellar systems which don't result in additional star 
formation) to influence the evolution of early-type galaxies (van 
Dokkum et al. 2003; Befl et al. 2003; van Dokkum 2005). How- 
ever, the dynamically hot environments in rich clusters which 
are the subject of our study are deleterious to the formation and 
survival of cold, bound-pairs of early-type galaxies - unless 
these systems arrive in the cluster as existing bound entities. 
It is not clear, therefore, that dry mergers can provide an ef- 
fective route to substantially increase the number of luminous, 
early-type SO galaxies within clusters. 

Unfortunately, the alternative mechanism for enhancing the 
luminosity of the bulge component - a starburst - also has 
strong observational evidence stacked against it. Surveys of 
star-forming galaxies in clusters using optical or UV star for- 
mation indicators have failed to detect galaxies with strongly 
enhanced star-formation which would have to exist to explain 
the growth of the bulge components of early-type galaxies in 
clusters at z ;S 0.5-1 (Balogh et al. 1999; Poggianti et al. 1999; 
Gerken et al. 2004). However, there is growing body of evi- 
dence that at least some of the galaxies in distant clusters may 
be undergoing bursts of star-formation, albeit ones which are 
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Fig. 1 .— The (B-R)-R color-magnitude diagram for (left) CI 0024+16 and (right) MS 0451-03. We identify 24-ij.m sources brighter than 200//Jy 
and compare these to the distribution for the optically-selected populations in these fields (for clarity, this has been represented by a smoothed 
density plot). Also shown are histograms showing slices through the distributions, comparing the number counts and colors for the mid-infrared 

sources to the optically-selected population. To emphasise the location of the sequence of early-type galaxies in the clusters, the (B-R) histogram 
is limited to R < 22. The mid-infrared selected population lies between the red and blue galaxy peaks in these fields, most likely because of the 
influence of dust on intrinsically blue, star-forming galaxies. For the mid-infrared sources in the color-magnitude diagram we have removed those 
objects morphologically classified as stars. All colors are measured in 2" diameter apertures, and we use the best_mag estimate of the total 
galaxy magnitude. 



heavily shrouded in dust. Small et al. (1999) used a deep VLA 
1 .4-GHz radio map to study a small sample of active galaxies 
within the core of die cluster C10939H-4713 (z = 0.41). Com- 
bining the radio data with near-infrared and optical morpho- 
logical information from the Hubble Space Telescope (HST) 
and ground-based spectroscopy, they found that over half the 
radio-emitting population in the core are dusty late-type galax- 
ies, presumably undergoing vigorous star formation. However, 
the spectral classification of these spirals placed them in the 
post-starburst class, and indeed all the post-starburst galaxies 
in this small region are radio-emitters. Dust has also been used 
to explain the unusual spectral properties of another class of 
galaxies found in distant clusters and the field: e(a) galaxies, 
which show enhanced Balmer absorption compared to normal 
star-forming galaxies (Poggianti et al. 1999). Poggianti & Wu 
(2000) and Poggianti, Bressan & Franceschini (2001) discuss 
models for these galaxies invoking age-dependent dust obscu- 
ration of the younger stellar populations - enabling significant 
activity to be hidden from view in these systems. If the passive 
lenticular galaxies found in local clusters, but absent from the 
equivalent rich environments at higher redshift, are the result of 
infaUing late-type galaxies undergoing dusty-starburst in high- 
z clusters, then a possible signature would be evolution in the 
total level of obscured star-formation in clusters out to z 1 . 

In principle, mid- and especially far-infrared/submillimeter 
observations give us a direct probe of the level of obscured ac- 
tivity in distant clusters. In particular, mid-infrared observa- 
tions with Infrared Space Observatory (ISO), and more recently 
Spitzer Space Telescope (SST), provide sensitive imaging capa- 
biUties which can trace dusty star formation in clusters out to 
z ~ 1 and beyond. Metcalfe, Fadda & Biviano (2005) sum- 



marise the results from ISO surveys of distant clusters, which 
have yielded a total of just ~ 40 cluster galaxies detected at 
15 iim across seven clusters between z = 0.18-0.56 (Due et al. 
2000, 2004; Fadda et al. 2000; Metcalfe et al. 2003; Coia et 
al. 2005a,b; Biviano et al. 2004). All these studies suggest that 
there is an increased level of mid-infrared activity in distant 
clusters, at levels above that suggested by UV/optical tracers 
of star formation. Submillimeter observations of more distant 
clusters have also hinted at possible enhanced activity in these 
environments (Best 2000; Webb et al. 2005). However, the in- 
homogeneous mix of coverage and depth in the samples cou- 
pled with the modest numbers of sources detected in any indi- 
vidual cluster mean that it has proved difficult to use these data 
to provide quantitative constraints on the origin and evolution 
of dust-obscured activity in distant clusters. 

The ^^'T's sensitive mid-infrared imaging capabilities pro- 
vide an unique opportunity to undertake complete and repre- 
sentative surveys of the obscured, active populations in distant 
clusters. To search for a population of mid-infrared sources in 
rich clusters environments, we have therefore used the Multi- 
band Imaging Photometer for Spitzer (MIPS) to detect 24-/xm 
emission from galaxies in two clusters at z ~ 0.5 covering a very 
wide range in environment from ~ 1 Mpc out to the turn-around 
radius (^ 5 Mpc) where the clusters merge into the surrounding 
field. These observations will provide measures of the level 
of obscured star-formation in these clusters, and so allow us 
to build up a reliable picture of the evolution of dust-obscured 
activity in clusters over the past 5 Gyrs. 

This paper presents a statistical analysis of the 24/im popu- 
lations in two z ~ 0.5 clusters. A subsequent paper (Geach et 
al. in prep) will discuss the properties of these sources in more 
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Fig. 2. — We show (R-K) versus (B-R) colors for 24^m-detected sources compared to the optical distribution in the CI 0024+ 16 (left) and 
MS 0451-03 (right) fields. The distributions are limited atR <24 and K <20 and we represent the optical sample's number density as a slightly 
smoothed grayscale. The dashed tracks trace the expected colors of galaxies across z = 0-2, with star formation histories similar to that of present- 
day Sdm~E galaxies, and the solid line and points show the location in color-color space of the range of spectral types at the cluster redshift. 
The arrow indicates the translation in color-color space corresponding to an increase in reddening of Av = 1. To perform a rough photometric 
cut to select 24^m sources at the cluster redshift, we define a region around the predicted colors of cluster members indicated by the rectangular 
selection boxes in both panels. To avoid contamination by infrared emission from AGN, we designate a sub-region which should contain the 
passive galaxies (dashed-line box) - 24/im sources in this region are not included in our analysis of the obscured star-forming populations. We 
show the efficacy of this selection by plotting the colors of known spectroscopically-confirmed 24/im cluster members, 85% of these fall within 
the selection box. 



detail using the available spectroscopic and morphological sur- 
veys of the clusters (Moran et al. 2006). The paper is organised 
as follows: we describe our observations and their reduction in 
§2, analyse these in §3 and discuss our results and present our 
conclusions in §4 and §5, respectively. Throughout, we adopt a 
geometry with 0^ = 0.3, ft a = 0.7 and Hq = 70 km s"' Mpc"' . 

2. OBSERVATIONS & REDUCTION 

The two clusters chosen for this study are unique in having 
panoramic HST imaging covering 25'-diameter fields - ex- 
tending from the cores out to the turn-around radii of the clus- 
ters (Treu et al. 2003; Kneib et al. 2003; Moran et al. 2005). 
These data have been used for weak-lensing analysis of these 
clusters, yielding 2-D maps of the dark matter distributions on 
~ 5-Mpc scales in the structures (Kneib et al. 2003). Panoramic 
studies of the galaxy populations in these clusters also bene- 
fit from extensive deep, ground-based optical and near-infrared 
imaging and spectroscopy (Moran et al. 2006 in prep). Al- 
though both clusters are relatively rich, they differ in their X-ray 
luminosities: CI 0024H-16 (z = 0.39) has a relatively modest X- 
ray luminosity, Lx 3.2 x lO'*^ ergs s"' (Treu et al. 2003), while 
MS 0451-03 (z = 0.55) is some 8x more luminous. This dis- 
tinction may lead to differences in the effectiveness of the var- 
ious processes influencing star formation (Treu et al. 2003), as 
traced by the distribution of the mid-infrared population. This 
will be useful in our subsequent detailed study to disentangle 
the potential mechanisms for triggering and suppressing star 
formation. 



2.1. Mid-infrared observations 

MIPS 24fim observations of the fields of C10024H-16 
(z = 0.39) and MS 0451-03 (z = 0.55) were obtained with 
SST in fixed-cluster offset mode on 2004 December 24-25. 
The observations of C10024H-16 are centered on 0026 35.70, 
+170945 (J2000); while those for MS 0451-03 are centered 
on: 045410.80, -030057 (J2000). The observations avoided 
the central ~ 5' x 5' of each cluster, which are part of Guaran- 
teed Time Observations (SST Program #83). For MS 0451-03 
we have acquired these data from the archive, and incorporated 
them into this work to provide coverage of the central region 
of the cluster The observations of the core of CI 0024+16 are 
not available at the time of writing, however there are existing 
mid-infrared observations from ISOCAM at 15/im (Coia et al. 
2005), which traces similar emission to that sampled by MIPS 
(namely emission from polyaromatic hydrocarbons). These al- 
low us to place limits on the mid-infrared population in the cen- 
tral regions of CI 0024+16. 

Each field was covered in a 5 x 5 grid (omitting the central 
pointing) yielding a 24-point mosaic. We adopted two 10-s cy- 
cles for our exposures producing a per-AOR pixel exposure of 
3 12.5s, and the total elapsed time for the AOR was 9.78 ks. The 
template AOR was repeated three times to provide redundancy 
and for the identification of cosmic rays and asteroids (the lat- 
ter can be common in 24/im MIPS data). The total exposure 
time per pixel is 938 s, which should achieve our goal of a 5-a 
detection limit of 200/zJy at 24/xm. 
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The basic-calibrated 24/im data provided by the Spitzer Sci- 
ence Center were first corrected for gradients, bright latents and 
the "jailbar" effect using calibration frames generated from the 
data. These were then mosaiked using the mopex package.^ 
Due to the 'first frame effect' (frames with a shorter exposure 
and a depression in response of up to 15%), the initial frame 
in each pointing was discarded, which improved flatness in the 
final mosaic. We implemented the the area-overlap interpola- 
tion method in the mosaicking procedure (Makovoz & Khan 
2005). The mosaics cover total areas of ~0.21 sq. degrees in 
both clusters. 

Radius (Mpc) 
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Fig. 3.— The lower panel shows the cumulative number of probable 
24/^m cluster members (brighter than 200/iJy) as a function of radial 
separation from their respective cluster centers. Probably cluster mem- 
bers are isolated based on their {R - K)-{B - R) colors (see Fig. 2). 
We correct the distributions for residual field contamination using a 
similar color-cut on the SWIRE ELAIS Nl survey (Lonsdale et al. 
2003, 2004, Surace et al. 2004), giving an estimated field correction 
of ~ 0.16 ± 0.02 arcmin"^. The estimated numbers of members could 
also be increased by ~ 18% to account for the incompleteness in our 
color-selection, although we have not done so here. We also plot the 
cumulative number of Ha-emitting (EW(Ha) >40A) cluster members 
in CI 0024+16 from the survey by Kodama et al. (2004), which exhibits 
a similar rise to that exhibited by the 24/im population in this cluster. In 
contrast, MS 0451-03 shows little evidence for a strong mid-infrared 
population. The top panel illustrates the variation in survey area with 
radius in our two fields, this shows that the slower increase in the cu- 
mulative 24/xm counts beyond 5 Mpc in CI 0024+16 is purely due to 
the decline in coverage in the outskirts of our survey. The errors on our 
observations are the combination of counting statistics and the mea- 
sured variance in the field number counts. The physical scale assumes 
an average redshift of z = 0.47. 

Source extraction was performed using SExtractor ver- 
sion 3.1 (Bertin & Amouts 1996) with the criteria that a source 
consists of at least 3 contiguous pixels (each pixel is 2.5" 
square) at > 2<t above the background. We measure 16-" diam- 
eter aperture fluxes, corresponding to ~ 3 x FWHM of the PSF 
(FWHM of ~ 5" at 24/im). Using a curve-of-growth analysis 
on bright isolated point sources, and our completeness simu- 
lations, we find that 16" apertures recover ^75% of the total 
flux and we therefore correct the resulting fluxes by a factor of 

^MOPEX is maintained by the Spitzer Science Center: http://ssc.spitzer.caltech.( 



1.33x to yield total fluxes. We visually inspect the extracted 
objects overplotted on the images, and remove sources which 
appear to be false, e.g. in the slightly noisier regions on the 
edge of the image. 

The flux detection limits and completeness of each mosaic 
was determined by inserting 10 sets of 10 simulated sources 
into the mosaics for a range of flux and determining their de- 
tection rate and recovered fluxes using identical extraction pa- 
rameters. The sources used in the simulations were high signal- 
to-noise composites formed by stacking a large number of real, 
isolated faint sources. From these simulations we estimate that 
the 80% completeness limits of our observations are 180/iJy in 
CI 0024+16 and 200AtJy in MS 0451-03. 

The archival mid-infrared data for the central region of 
MS 0451-03 (Program #83), were obtained from the Spitzer 
archive. The observations were taken on 2004 September 23, 
centered on 045410.8, -03 0057 (J2000). We use die post- 
basic-calibrated data, and subject it to the same extraction crite- 
ria as described in §2.1. These data are slightly shallower than 
our mosaic, and we take this into account in our subsequent 
analysis. 

2.2. Archival optical & near-infrared imaging 

We use existing deep optical and near-infrared imaging to 
obtain colors for the 24/xm sources in these regions. For both 
clusters we use Subaru SuprimeCam B- and /?-band observa- 
tions taken for the PISCES survey (Kodama et al. 2005). The 
CI 0024+ 16 field was observed on 2002 September 6 under 
good conditions. The seeing was ^ 0.7-1" for the R- and ^ 1- 
1 3" for the £-band, with exposure times of 5,280 s and 3,600 s 
respectively, reaching a depth of /? ~ 27 (we use Vega-based 
magnitudes throughout). The observations and data reduction 
technique are described in Kodama et al. (2004). MS 0451-03 
was observed on 2001 January 21-22, again in good conditions 
wifli 1.0" seeing in tiie fi-band and 0.8" in flie /?-band. Total 
exposure times were 7200 s in B and 4800 s in R, again yielding 
photometry for objects as faint as ^ 27. 

Panoramic near-infrared /T-band imaging of both clusters is 
also available with WIRC (Wilson et al. 2003) on the Palo- 
mar Hale 5.1-m telescope. These data comprise a 3 x 3 mo- 
saic of WIRC pointings, providing a contiguous observed are of 
26' X 26' centered on each cluster Full details of the observa- 
tions and data reduction are published elsewhere (CI 0024+16: 
Smith et al. 2005; MS 0451-03: Snuth et al. 2006, in prep). 
Point sources in the final reduced mosaics have a fuU width half 
maximum of 0.9" and 1 .0" in CI 0024+16 and MS 0451-03 re- 
spectively. The data reach a 5-(t point sources detection thresh- 
old of = 19.5 and = 20.0 in CI 0024+16 and MS 045 1-03 
respectively. The nominal 24/tm coverage of the two clusters is 
-0.21 degrees^. In the case of MS 0451-03, the archival GTO 
data contribute an extra ~29 arcmin^. When the optical and 
near-infrared coverage is taken into account, the total coinci- 
dent coverage of 24/im, B-, R-, A'-bands is 0.165 degrees^ and 
0.184 degrees^ in CI 0024+16 and MS 0451-03 respectively. 

We detect objects and extract photometry from the optical 
frames by running SEXTRACTOR in two-frame mode such that 
detections were made in the R-band image and measurements 
made at identical pixel locations in R and B (the frames were 
previously aligned to good accuracy). We catalog all sources 
with 3 contiguous pixels (the pixel scale is 0.204") at least 2-a 
above the background, and lay down 2" diameter apertures to 
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measure colors in the BRK-hand frames. 

3. ANALYSIS & RESULTS 

The MIPS observations of our two clusters yield 986 sources 
in the CI 0024+16 field and 1071 in MS 0451-03, both with 
24/xm flux densities above S24fj.m > 200y^Jy (unless otherwise 
stated, quoted fluxes are always the corrected 16" aperture val- 
ues). In the case of MS 0451-03 the count includes sources in 
the archival region in the core. We assess the false detection 
rate by running our algorithm on the inverse of the data, and 
detect ~ 20 sources in each frame, all with fluxes in the range 
200-330/iJy, representing ~2% of the detected sample. Our 
analysis requires matched optical R <24 identifications and so 
we expect the false rate will be much less than 1%. As we 
demonstrate below, the achieved source surface density is be- 
low 1 object per 40 beams, the classical definition of a con- 
fused map (Hogg 2001), and so our maps are not expected to 
be confused at this depth. 

We next use the deep optical and near-infrared imaging of 
these fields to investigate the photometric properties of mid- 
infrared sources. The influence of false detections described 
above is minimised in the main part of our analysis, since we re- 
quire the mid-infrared sources to have optical and near-infrared 
counterparts. 

3.1. Optical properties 

Positional matching of multi-wavelength data sets has been 
a long standing problem in astronomy, and can be particularly 
troublesome when there is a large disparity in the resolution 
and sampling in two datasets. We find that a simple nearest- 
neighbor match is not adequate to pair mid-infrared sources 
with optical counterparts - the method can fail to match com- 
plicated interacting systems for example. Instead, we apply 
the technique of de Ruiter, Arp & Willis (1977) who use a 
Bayesian estimator for the probability, p(id | r), that a nearby 
source is a true match and not a chance unrelated object. In 
Appendix A, we briefly outline the key elements of the method, 
but refer the reader to de Ruiter, Arp & Willis (1977) for a thor- 
ough derivation. The result of our matching analysis is a Ust of 
probable optical counterparts for the 24 /xm sources in our cat- 
alogs. We identify 611 and 650 counterparts to 24/im sources 
which are brighter than R = 24 and 200 /uJy in CI 0024+16 and 
MS 045 1-03 respectively. In MS 045 1-03 this includes sources 
in the archival region. 

In Figure 1 we plot the (B-R) versus R color-magnitude di- 
agram for each cluster, comparing the distribution of colors for 
the mid-infrared sources and their apparent magnitudes with 
the optical population in these regions. The first thing to note is 
the broad similarity in the distributions on the color-magnitude 
plane of the 24^m sources in the two fields. Looking at the 
color distributions in more detail, the median (B-R) color of 
mid-infrared sources at /?< 22 is 1.75 and 1.80 for CI 0024+16 
and MS 0451-03 respectively. The 24-/im counterparts exhibit 
a broad peak in (B-R) color, which lies between the peaks 
of red and blue galaxies (Fig. 1). This association of 24/Ltm 
sources with galaxies having transition colors is intriguing. It 
may simply reflect the fact that these galaxies are dustier ex- 
amples of the general blue star-forming field population, where 
dust reddening produces somewhat redder colors. Or it may be 
evidence that some of these galaxies are part of an evolutionary 
sequence connecting the blue star-forming population, and the 
passive types inhabiting the red-sequence. 

In Figure 2 we plot the (R-K) versus (B-R) colors of 



the mid-infrared sources in CI 0024+16 and MS 0451-03 com- 
pared to the optically-selected populations in the two fields. 
We also indicate the expected colors of cluster members with 
a range of spectral types (King & ElUs 1985). These model 
colors provide a rough guide to the relative level of current to 
past star-formation in the galaxies and the reader should note 
that the influence of dust extinction will move galaxies parallel 
to this sequence, further complicating any interpretation of the 
star-formation histories of these galaxies. 

Looking at the panel for CI 0024+16 in Figure 2, we see a 
clear ridge of 24-^m sources with colors comparable to those 
expected for cluster galaxies. To confirm this association, we 
indicate on the plots the locations of spectroscopicaUy con- 
firmed 24^m cluster members from the surveys of Czoske et 
al. (2001) and Moran et al. (2006), and our on-going spectro- 
scopic follow-up of our 24/xm sample (Geach et al. in prep). 
There are 45 confirmed 24/im members in CI 0024+16 and 7 
in MS 0451-03 and the vast majority of these lie close to the 
predicted color-color locus for cluster members, and have col- 
ors consistent with star-forming galaxies. However a few fall 
in the region populated by passive, early-type cluster members, 
as expected due to the likely presence of active galactic nuclei 
in some of these galaxies. Nevertheless, it is clear that the in- 
formation from this {B-R)-{R-K) color-color plot can be used 
to crudely isolate likely 24-/um cluster members to statistically 
gauge the size of this population in the two clusters. 

The 'transitional' colors of the majority of 24/tm sources 
hinted at in the color-histograms in Figure 1, is also reflected 
in the distribution of points around the spectral classes in Fig- 
ure 2. The colors of the bluer, probable 24/im cluster members 
suggest these are moderately star-fornoing galaxies (matching 
the continuum shapes of present-day Sbc/Scd galaxies). How- 
ever, these could be more actively star-forming systems with 
additional reddening, ~ 0-1, and there are a small number 
of galaxies with very red colors which suggest very significant 
reddening. Ay ~ 2 mags. To select the majority of infrared clus- 
ter members, we apply a photometric cut in color-color space 
around the sequence of known cluster members (Figure 2) and 
then account for the residual field contamination using a simi- 
lar color-selection applied to a blank field 24/tm sample. Based 
on the existing spectroscopic samples of 24/im sources in the 
two clusters, we estimate that our color-selection criteria iden- 
tify 39/45 and 5/7 members in CI 0024+16 and MS 0451-03 
respectively, or an average completeness of 85%. We note that 
the fact that the reddening vector runs approximately paraUel 
to our selection box in Figure 2 means that it should contain a 
significant fraction of the more reddened sources. 

In addition, in this work we are investigating the nature of the 
star-forming cluster population, so we must be aware of con- 
tamination by infrared activity caused by AGN in passive (and 
also star-forming) galaxies. Therefore we define a sub-region 
in the color-space around the likely color of cluster E/SOs (Fig- 
ure 2), and do not include the galaxies within this region in our 
analysis of the obscured star-forming population. 

The 24/im galaxy population in the MS 0451-03 field does 
not exhibit such a pronounced ridge with colors matching those 
expected for cluster members, as CI 0024+16. Nevertheless, 
there is a weak overdensity, which is most obvious in a red 
clump associated with passive, early-type cluster members. 
This clump is around 0.2 mags redder in (B-R) and (R-K) 
than the equivalent feature in CI 0024+ 1 6, as expected from the 
^-correction for passive galaxies between z = 0.39 and z = 0.55. 
Again, the small number of spectroscopicaUy confirmed 24/tm 
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members fall close to the locus expected for the colors of clus- 
ter members and we therefore use this color-color selection in 
the subsequent sections to investigate the properties of this pop- 
ulation. 
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Fig. 4.— The infrared luminosity function for CI 0024+16, based on the 
photometrically-selected cluster members and plotted with and without 
correction for residual field contamination. The total infrared luminosi- 
ties are derived by extrapolating the monochromatic MIPS 24/xm flux, 
assuming the local 15^m-far-infrared correlation holds at z ~ 0.5 (see 
§3.3). The estimated number of 24/im cluster members in MS 045 1-03 
is insufficient to provide a similar plot for that cluster. For compari- 
son, we plot the luminosity function of the 60;um IRAS Point Source 
Catalog Redshift (PSCz) survey from Takeuchi et al. (2003), with an 
arbitrary normalisation. The cluster luminosity function is similar in 
shape with that of3 local field galaxies, given the large uncertainties. 

3.2. Counts of mid-infrared sources 

To quantify the number of 24/im sources associated with 
cluster members we first compare the 24yum number counts for 
each cluster to counts from the literature for large area 24-/zm 
field surveys (e.g. Marleau et al. 2004; Papovich et al. 2004; 
Chary et al. 2004). The counts in the literature all agree well 
in the flux density regime covered by our data and we use the 
model presented in Marleau et al. (2004) to estimate the num- 
ber of field 24/im galaxies expected with 524^m = 200-1000/iJy 
in our fields: dN/dS = (416 ± S)^-^ ''^" ', which is valid for the 
range 200-900/iJy, we extrapolate the bright end to lOOO/xJy. 
For this comparison we also correct our counts for complete- 
ness in the range 200-300/iJy (our simulations show that our 
observations are ~100% complete at 300/zIy). 

The predicted numbers of 24 fim sources to our flux limits 
over flie CI 0024H-16 and MS 045 1-03 fields are 1280±160 and 
939±120 respectively, down to 200/Jy. These should be com- 
pared to the observed numbers sources of 1032 and 935, re- 
spectively. Based on this comparison we can see that the clus- 
ter populations do not represent a significant excess over the 
large field population across the fuU survey area and down to the 
flux limits of our survey. This is not particularly surprising as 
the survey fields are large and our sensitivity limit is sufficient 
to detect sources in a large volume out to high-redshifts. In- 
stead, if we restrict the comparison to the central regions of the 
clusters, where the contrast of the cluster over the background 
might be largest, we see that at projected radii of < 2 Mpc there 



is a slight excess of mid-infrared sources compared to the field 
counts in CI 0024H-16: 134 ± 12 observed sources compared to 
115 ± 14 predicted for the field, suggesting a modest cluster 
population. While in MS 0451-03, including the archival data 
for the core, and we find that within 2 Mpc there are 104 ± 10 
sources compared to 108 ± 13 predicted by the model - show- 
ing no detectable overdensity in this cluster even in the central 
regions. 

The lack of a detectable excess in the raw counts of 24/xm 
sources in the cluster fields does not rule out a significant clus- 
ter population given the large uncertainties in the contribution 
from field sources. To provide a more sensitive test we can ex- 
ploit the optical-near-infrared photometry to remove the bulk 
of the background field contamination and so reduce the uncer- 
tainty in the field correction. In Figure 3 we plot the cumu- 
lative number of mid-infrared galaxies as a function of clus- 
tocentric radius for those galaxies with optical-near-infrared 
counterparts whose colors Ue close to the ridge-lines of clus- 
ter members seen in Figure 2 (restricted to those with R < 24). 
These are corrected for residual field contamination by using 
the optical/near-infrared photometry of mid-infrared sources 
from the SWIRE ELAIS Nl survey (Lonsdale et al. 2003, 2004; 
Surace et al. 2004). To achieve this we translate our (B-R)- 
(R - K) selection box in Figure 2 to the equivalent region in 
{r-z)-{g-r) space for the SWIRE sample using model spec- 
tral energy distributions (Fig. 2) to transform between the var- 
ious colors. Extrapolating the number counts of these color- 
selected sources to our 200/iJy limit (by fitting a power-law, 
dN/dS oc S~", to the counts) yields an average surface den- 
sity of mid-infrared field sources with /? < 24 of 0.16 ±0.02 
arcmin"^ and we correct our counts using this. We also make 
use of the very large SWIRE survey area to estimate the fluc- 
tuations in the field correction in the area of each of our radial 
annuli in Figure 3, by determining the fluctuation in the field 
counts as a function of survey area. This variance is included in 
the uncertainty in the field correction used in Figure 3. 

By selecting 24/im sources with optical-near-infrared colors 
similar to those expected for cluster members, we have sub- 
stantially reduced the field contamination in our sample. Fig- 
ure 3 demonstrates a significant 24/xm cluster population in 
CI 0024+16, with 155±18 sources across our 25' x 25' field 
(where the uncertainty includes the effects of clustering on the 
field correction). As can be seen the excess population is dis- 
tributed across the whole field, with the numbers continuing to 
increase to ~ 5 Mpc where the coverage of our data begins to 
decline. In contrast, in MS 0451-03 we find a much weaker 
excess, 28 ±17, which is only marginally significant. Never- 
theless, we note that our existing spectroscopic survey of this 
cluster (Fig. 2) confirms that there are 24-/xm-detected cluster 
members, but the population appears much less numerous than 
in CI 0024+16. In both cases, the estimates of the total clus- 
ter populations could be increased by ~ 18% to account for the 
incompleteness of our color-selection, although to be conserva- 
tive we choose not to apply this correction. 

While it is clear that there are populations of 24/tm clus- 
ter members, we also need to consider lensing of background 
sources by these massive clusters, which may affect the pre- 
dicted numbers of cluster members. Lensing has the effect of 
a boosting of flux and dilating the projected area. To demon- 
strate this we model the clusters as simple isothermal spheres. 
In this case a background source will suffer an apparent de- 
flection in radial direction by a ~ 20"(Tiooo> where ctiooo is 
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the cluster velocity dispersion in units of 1000 km s '. We 
adopt a = 1 150km s-i for CI 0024+16; a = 1290km s"' for 
MS 0451-03 (Trentham et al. 1998). At a radius r, the source- 
plane area of an annular bin of width dr will be 27r(r- a)dr, and 
the flux of a lensed source will be boosted by a factor r/{r-a). 
We integrate the lensing prediction for our color-selected sam- 
ple (taking a = 1200 km s~' - though the actual values may be 
lower, so this should be a conservative estimate) as a func- 
tion of radius, finding that it could contribute < 20 sources to 
our samples, hence in CI 0024+16 the excess signal we see is 
greater than can be explained by lensing. The situation with 
MS 0451-03 is more problematic, although again we note that 
the spectroscopic identification of a handful of 24/im cluster 
members confirms that, as with CI 0024+ 16, there is a popula- 
tion of mid-infrared sources associated with the cluster 
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Fig. 5.— a comparison of the star-formation rates derived from Ha 
(Kodama et al. 2004) and total-infrared luminosities for individual Ha- 
detected cluster members in CI 0024+16. The ratio SFR[IR]/SFR[Ha] 
provides an estimate of the level of obscuration of the starforming re- 
gions within these galaxies. We see a typical offset of a factor of 5 x 
between the infrared and Ha estimates (the latter have not been cor- 
rected for extinction), although there is a wide-scatter. We also plot 
data from the literature of similar studies of infrared-luminous field 
galaxies from Flores et al. (2004), Franceshini et al. (2003) and Dop- 
tia et al. (2002) corrected to our cosmology. It appears that the level of 
obscuration of star-forming infrared cluster galaxies (in CI 0024+16) is 
similar to that of mid-infrared selected field galaxies at this epoch. We 
do not plot individual error bars for clarity, but show a representative 
error bar. Finally, to probe obscuration from faint mid-infrared sources, 
we stack the 24/im data for Ha emitters which are not detected in the 
MIPS image yielding a 5a statistical detection. This confirms the gen- 
eral trend seen above our nominal luminosity limit that there appears 
to be an increase in the level of obscuration (SFR[IR]/SFR[Ha]) with 
infrared luminosity. 



3.3. Star formation 

To calculate the star-formation rates (SFRs) for our mid- 
infrared detected galaxies, we begin by estimating the total in- 
frared luminosity over the restframe wavelength range 8-1000 
fim. Since we only have MIPS 24/im photometry, we have to 
determine the correction factor required to convert our 24/im 
luminosities to total infrared luminosities. To achieve this we 
assume the correlation between 15/im (corresponding to ob- 
served ~ 24 jim at z ~ 0.5) and total-infrared luminosity seen 
for local infrared galaxies (Chary & Elbaz 2001) continues to 
hold at z ~ 0.5. We then follow a similar technique to that of 
Bell et al. (2005) and use the spectral energy distribution (SED) 
templates from Dale & Helou (2002) to estimate the ratio of 
the observed-frame 24-/tm luminosity to the 8-1000/im lumi- 
nosity. We calculate the ratio of total-infrared to observed 24- 
/tm luminosities for each SED template, taking the mean value 
of all the ratios as our correction factor, and the range from 
the most active to the most quiescent SEDs as a conservative 
estimate of the systematic uncertainty. The correction factors 
are 16 ± 2 and 16 ± 3 for CI 0024+16 and MS 045 1-03 respec- 
tively. We compare our conversion factors to the relation found 
by Chary & Elbaz (2001) by calculating the mid-infrared to to- 
tal infrared ratio with the ISOCAM UN 3 (15/im) filter for the 
Dale & Helou templates at z = 0.1. The Chary & Elbaz conver- 
sion for local LIRGs corresponds to Lir ~ 1 1^4 x while 
our estimate gives a prefactor of 7 ± 3 over the full range of 
LIRG templates. Thus these two calibrations are consistent, al- 
though the reader should note our estimates are roughly 40% 
fainter than would be given by Chary & Elbaz (2001). 

Our 5 G 24-/im flux limit of 200/iJy corresponds to average 
total infrai-ed luminosities of 6 x lO'^L© and 12 x lO'^L© in 
CI 0024+16 and MS 0451-03 respectively, hence the bulk of 
the populations we are detecting are Luminous Infrared galax- 
ies (LIRGs) with L/« > 10" L©. These luminosities translate 
into SFRs of ^ 10 M© yr"' and ~ 20 M© yr~' assuming the far- 
infrared star-formation calibration given by Kennicutt (1998) 
with a Salpeter IMF and a mass range of 0. 1-IOOM©. 

In Figure 4 we plot the mid-infrared luminosity function for 
the CI 0024+16, based on the color-selected sample out to a 
radius of 5 Mpc. To field correct the luminosity function in 
CI 0024+16, we integrate the extrapolated power-law fit to the 
color-selected 24/im sources in the SWIRE ELAIS Nl field, 
from §3.2. We use the scatter in the normalisation of dN /dS 
in independent 25' x 25' regions within the SWIRE field, com- 
bined with the counting errors in each bin, to give the uncer- 
tainty in the field-correction. This plot demonstrates that the 
field-corrected luminosity function in CI 0024+16 is similar in 
form to that seen for 60/im-selected populations at low-redshift 
(Takeuchi et al. 2003). 

While there are too few sources in MS 0451-03 to allow us 
to make an equivalent plot to Figure 4 for that cluster, we can 
use the CI 0024+16 luminosity function to estimate how many 
sources we would have detected if we place the CI 0024+ 16 
population at z = 0.55. Applying a factor to scale for the effec- 
tive areas of our surveys in CI 0024+16 and MS 0451-03, we 
estimate that we would detect 69 ±21 galaxies above the lu- 
minosity hmit, 12 X 10'°Lo, of the MS 0451-03 survey. This 
compares to our estimate of 28 ± 17 from our MS 045 1-03 cat- 
alog, indicating a roughly 2-3 x difference in the numbers of 
mid-infrared sources in the two clusters. This suggests that, al- 
though the dearth of mid-infrared sources in MS 0451-03 may 
be partly due to sensitivity limitations, there may be a real dif- 
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ference in the mid-infrared populations in the two clusters. 

3.3.1. Mid-infrared versus Ha tracers of activity 

We can also compare the properties of our mid-infrared se- 
lected sample in CI 0024H-16 with the Ha survey of this cluster 
by Kodama et al. (2004). Kodama et al. obtained a deep im- 
age of the cluster in a narrow-band filter centered on redshifted 
Ha. Galaxies with excess emission in this narrow-band filter 
were identified as cluster members with Ha emission. They es- 
timate their sensitivity limit corresponds to a star formation rate 
of ^ 1.5 M0 yr"' - significantly below the estimated SFR limit 
of our mid-infrared survey, ~ 10 yr"'. 

To start with, we compare in Figure 3 the cumulative radial 
profile of the mid-infrared sources with the Ha-emitting galax- 
ies from Kodama et al. (2004). As can be seen the two trac- 
ers show very similar trends, with a steep rise from 1 Mpc out 
to ~5 Mpc, beyond which our coverage declines. The reader 
should also note that the lack of mid-infrared coverage for the 
core of the cluster will contribute to the steeper drop of the 
24/im counts in the innermost bin in Figure 3. Surprisingly, 
Figure 3 also shows there are similar numbers of Ha emitters 
and 24/Ltm-detected members, even though there is a signifi- 
cant difference in the relative sensitivity of the two surveys. We 
therefore next compare the estimates of the SFR for individ- 
ual galaxies from Ha and the mid-infrared, to relate these two 
star-formation indicators for confirmed cluster members. 

In Figure 5 we compare the star-formation rates for Ha 
narrow-band- selected sources and their 24/xm counterparts^. 
From this figure, it can be seen that the star formation rates 
for the mid-infrared-detected members, derived from the Ha 
emission line (but not corrected for extinction), appears to un- 
derestimate the SFR from the mid-infrared by factors ^ 5 x . In- 
cluding the standard assumption of 1 magnitude of extinction 
and correcting for a 30% contribution to the measured fluxes 
from [Nil] (see Kodama et al. 2004) would reduce this offset to 
a factor of 2.5 x . However, such a simplistic correction would 
miss the fact that the discrepancy in SFRs appears to increases 
with infrared luminosity, such that the most active systems are 
also the most obscured. 

To test for biases due to our 24/xm sensitivity limit, we stack 
thumbnail images from the MIPS image of CI 0024H-16 for Ha 
emitters from Kodama et al. (2004) which are not individ- 
ually detected, i.e. 524^in < 200 /xJy. We find a ~5-cr de- 
tection with >S'24;nm = 57 ± 11 fxJy. In terms of star-formation, 
this corresponds to '~3 Mq yr~^ We plot this point on Fig. 5, 
with the corresponding average Ha derived SFR for the stack, 
^ 1 Mq yr"' . This point seems to be roughly consistent with the 
trend seen, that obscuration (i.e. SFR[IR]/SFR[Ha]) increases 
with the level of infrared activity. 

We estimate that the correction required to bring the Ha de- 
rived SFRs in line with the infrared SFRs corresponds to a red- 
dening of Ay ^ 2.5, slightly larger than the value used to ex- 
plain the broadband colors of these galaxies on the (R-K)- 
(B-R) color-plane. This reddening estimate will also be af- 
fected by the choice of adopted SED when deriving total in- 
frared luminosities, although in this work we have chosen a 
conversion which represents infrared galaxies of a wide range 
of activities, and so our estimates are conservative. We also note 
that there is a factor of 40 range in infrared to Ha derived star 
formation rates within our sample which is larger than expected 
from our measurement errors and indicates that there may be a 

'The optical star-formation rates are derived from the Ha and adjacent [N u 



wide range of extinction amongst the individual LlRGs - in part 
probably due to the patchiness of the dust within these galaxies. 

We also compare our relation between mid-infrared- and Ha- 
derived SFRs to that seen in samples of field galaxies from the 
literature: a sample of LlRGs in the range 0.2 < z < 0.7 from 
Flores et al. (2004), 15/L(m detected ISOCAM sources from the 
Hubble Deep Field South at (0.2 < z < 1-5) (Franceshini et al. 
2003), and IRAS 60/im warm infrared galaxies from Doptia et 
al. (2002). The typical level of SFR[IR]/SFR[Ha] is similar to 
Hq/MIR detected objects in CI 0024H-16, suggesting that these 
cluster galaxies are similar in nature (at least in terms of their 
mix of star-formation and obscuration) to star-forming LlRGs 
in the field at similar epochs. 

3.4. Total cluster star-formation rates 

We now estimate a total star-formation rate for the clusters 
using the typical SFR of the likely star-forming mid-infrared 
members. Since we lack observations in the central ^ 1 Mpc in 
CI 0024h- 16 due to the unobserved GTO region, we also need to 
correct for this missing population using the 15/im ISO obser- 
vations from Coia et al. (2005). 

3.4.1. Cl0024-\-16 

As shown by the radial number counts presented in Figure 3, 
the 24/im members are widely distributed across the clusters 
out to a projected radius of ~2Mpc (corresponding to ~ /?200> 
Carlberg et al. 1997. Note that throughout, we use projected 
distances) and beyond. To compare the mid-infrared-estimates 
of the total star formation rates within representative regions 
of different clusters we therefore adopt a 2 Mpc radius for our 
calculations, which roughly matches the size of the surveyed 
regions in some of the earlier small-field surveys. 

We integrate over the luminosity function in Figure 4 for 
sources with projected radii < 2 Mpc, correcting for the ex- 
pected contamination from field galaxies using the SWIRE data 
for the ELAIS Nl field (see §3.2). To account for the miss- 
ing sources in the unobserved GTO core region we need to add 
the SFR in the 15 members in this region found by Coia et 
al. (2005), after converting their 15/xm fluxes to star-formation 
rates using an identical method to that presented in §3.3. 

Our derived total-infrared luminosities are lower than those 
derived by Coia et al. (2005) by a factor of ^2. This difference 
can be attributed to the different methods of conversion used. 
Coia et al.'s method involves extrapolating their observed 15/ttm 
photometry to 15/im in the rest-frame, using an SED fitted to 
their optical and mid-infrared photometry. This is then con- 
verted to a total infrared luminosity using the empirical conver- 
sion of Chary & Elbaz (2001) described above. In contrast our 
observations correspond to ~15/im in the rest-frame. Hence 
some of this offset is due to the differences in conversion from 
15/im luminosity to total infrared emission discussed in §3.3 
and the remainder is due to the choice of SED by Coia et al. 
(2005), which overpredicted the observed 24/Ltm fluxes. As we 
have stated earlier, our calculations should yield conservative 
limits on the SFR. We note that, with our total-infrared cali- 
bration, the 5a detection limit of the Coia et al. (2005) ISO 
observations corresponds to ~ 3 x lO'^L©. 

With our conversion, we find three galaxies in the Coia et al. 
(2005) sample that exceed our luminosity limit of 6 x lO'^L©, 
and so we can add their SFR to our total within 2 Mpc of 
the core. We thus find a total star-formation rate for the 40 

48,6583 emission lines after Kennicutt, Tamblyn & Condon (1994). 
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sources within 2Mpc to be 1000 ± 210 M0 yr~'. The error 
is derived by boot-strap resampUng the flux distribution of 
sources within 2 Mpc, and integrating the new luminosity func- 
tion. The median SFR per galaxy in our sample within 2 Mpc 
is ^ 16 M0 yr"'. We note that correcting for the estimated in- 
completeness in our color-selection of cluster members would 
increase our estimate of the total SFR by a factor 1.07, which 
corresponds to a change less than its quoted uncertainty. 

For comparison, we also note that the Ha observations of 
Kodama et al. (2004) detect 100 galaxies within a similar re- 
gion to that discussed here. They find a cumulative SFR of 
^ 47OM0yr"', above their approximate sensitivity limit of 
SFR of 1 .5 Mq yr"\ or about 4.7 M© yr"' per galaxy. Restrict- 
ing the Ha sample to the > 10 Mq yr"' limit of our mid-infrared 
survey, we find that the Ha survey would yield an integrated 
SFR of ~ 220 Mq yr~'. This confirms the effect seen in Fig- 
ure 4, that although Ha and mid-infrared surveys are detect- 
ing a similar population, the Ha tracer severely underestimates 
the underlying activity in the most active sources in the clus- 
ter: these starbursting systems are producing dust at a more 
copious rate and therefore optically obscured. In addition, we 
can state that the mid-infrared survey detects the bulk of the to- 
tal star-formation activity within this region, requiring only a 
~20% correction for the star-formation activity in sources with 
1-10 M0yr-'. 

3.4.2. MS 045 1-03 

We have seen that, compared to CI 0024H-16, this cluster ap- 
pears to be deficient in mid-infrared sources (at least down 
to our luminosity limit). Nevertheless, even a small excess 
of mid-infrared sources in the cluster could contribute a non- 
negligible star-formation rate to the overall activity of the clus- 
ter. We therefore estimate the total star-formation rate within 
~2Mpc using the small number of excess sources identified 
using our {R-K)-{B-R) color-color selection. Applying 
the same approach as used for C10024H-16, we find a total 
star-formation rate within a 2 Mpc radius of the cluster cen- 
ter of 200 ± 100 Mq yr"'. The median SFR within 2 Mpc is 
~ 35M0yr"' per galaxy (reflecting the brighter luminosity 
limit in this cluster) . 

In §3.3 we estimated the number of mid-infrared sources we 
would detect if MS 0451-03 had the same mid-infrared lumi- 
nosity function to C10024H-16, taking into account the slightly 
different areal coverages between the clusters. Similarly, we 
can estimate the total star-formation rate in MS 0451-03 if we 
reached a luminosity limit identical to C10024H-16. Assuming 
the faint end of the luminosity function in MS 045 1-03 fol- 
lows a similar shape to C10024-H16, we estimate that the to- 
tal star-formation rate down to 6 x lO'^L© in this cluster is 
< 460 M0 yr"'. Although the small number of excess sources 
in MS 0451-03 can lead to a relatively high star-formation rate 
(compared to what might be found using an optical tracer for 
example), there is a clear difference in the population of mid- 
infrared sources in the two clusters - but what is the physical 
cause of this? 

Although the clusters are of similar mass, and at a similar 
redshift, they differ strongly in their intra-cluster environments 
and dynamical status (Treu et al. 2003). CI 0024H-16 is dynam- 
ically active (Czoske et al. 2001; Kodama et al. 2004) and this 
may provide the impetus for the triggering of star-formation 
via mergers and interactions of gas-rich spirals within the clus- 
ter - even in the apparently high-density core. The intraclus- 



ter medium (ICM) of MS 0451-03 is much hotter and denser 
(by nearly an order of magnitude) than C10024H-16, and this 
will have an impact on the radii that processes such as star- 
vation and ram-pressure stripping operate (Treu et al. 2003). 
In MS 0451-03 the hot ICM will be much more effective at 
starving galaxies of their gas reservoirs at a larger clustocentric 
radii than in C10024H-16. The relative dearth of mid-infrared 
sources in MS 045 1-03 might then suggest that the active re- 
gions within the mid-infrared population are comparatively sen- 
sitive to these gas removal mechanisms. Unlike C10024H-16, 
unfortunately there are no published optical SFR studies of 
MS 045 1-03 to confirm a wholesale decline in the SFR in 
this cluster. Nevertheless, the small excess of mid-infrared 
sources in MS 0451-03 show that it hosts quite significant star- 
formation, so it is unclear at this stage exactly what physical 
processes control the distribution of star-forming galaxies in 
different types of clusters. We will address the issue in our next 
paper where we study the properties of the 24/um cluster mem- 
bers in more detail using the available spectroscopy, dynamics 
and HST imaging (Geach et al., in prep.; see also Moran et al., 
in prep.). 

4. OBSCURED STAR-FORMATION IN CLUSTERS OUT TO Z ~ 1 

At the present-day, virialised regions of the Universe are 
dominated by passive galaxies: ellipticals and lenticulars, with 
the main contribution to the global star-formation rate den- 
sity at z = coming from late-type spirals in low-density en- 
vironments. However, at higher redshifts the high-density envi- 
ronments within clusters show increasing numbers of actively 
star forming galaxies, possibly reflecting a similar increase in 
the frequency of activity to that seen in the surrounding field. 
Moreover, it is clear that there is significant hidden star forma- 
tion in clusters in this redshift range: observations by ISO and 
now our new Spitzer observations reveal populations of infrared 
galaxies in clusters at low and intermediate redshift, with SFRs 
much higher than would be measured using optical tracers such 
as Ha or [Oii]A3727, due to the extinction effect of dust. 

Given this new information, what is the evolution of the total 
star-formation rates in massive clusters - i.e. how does the in- 
crease in activity in high-density environments at high redshift 
behave in the mid-infrared? To compare the total mid-infrared 
derived star-formation rate in clusters from z = 0-1 we use the 
total SFR within ~2Mpc and normalize by the best-estimate 
mass of the cluster (for example, see Kodama et al. 2004). For 
C10024H-16 and MS 0451-03 these are 6.1 x lO'^'M© (Kneib et 
al. 2003) and 15 x lO'^^M© (La Roque et al. 2003) respectively, 
where the masses are within 2 Mpc of the clusters' centers. 

To build up a history of star-formation in clusters out to 
z ~ 1, we compare our results to previous ISO and IRAS stud- 
ies of the clusters: Perseus (A 426), A 1689, A 370, A 22 18 and 
J1888.8CL (CI 0054-27), summing the known cluster mem- 
ber's star formation rates down to a far infrared limit 6 x lO'^L© 
- the limit of the C10024H-16 observations. To ensure correct 
comparison of the rates, if necessary we re-derive the SFRs 
using our employed calibration from Kennicutt (1998), with a 
standard Salpeter IMF with a mass range O.l-IOOM©. For the 
ISO observations, we also re-derive the far-infrared 8-1000 /zm 
luminosity using the method we present above, using the band- 
pass of the ISO LW3 filter to convert mid to total-infrared. The 
ISO observations of these clusters are summarized in the re- 
view by Metcalfe et al. (2005). Note that the previous ISO 
and IRAS observations have concentrated on the core regions 
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of clusters (typically within 1 Mpc), so we restrict our star- 
formation integration to within a radius of 2 Mpc in CI 0024+16 
and MS 045 1-03, and compare to estimates of the total star- 
formation within the equivalent radius in the other clusters. In 
the case of MS 0451-03, we estimate an upper limit to the star- 
formation rate extrapolated to a luminosity limit of 6 x 10'" Lq 
by using the method described in §3.4.2, assuming the faint 
end of the luminosity function follows a similar shape to that 
in CI 0024+16. We estimate that the upper limit to the to- 
tal star-formation rate down to 6 x 10^" L© in MS 0451-03 is 
<460MQyr"'. 




Redshift 

Fig. 6.— The variation in the mass-normalised star-formation rates in 
clusters out to z ~ 0.5. The star formation rates are from the mid- 
infrared populations within ~2Mpc and these are normalised to the 
best estimate of the cluster mass, which is derived via lensing esti- 
mates, or from the X-ray luminosities. We also plot an evolution- 
ary model for the counts of star-forming ULIRGs from Cowie et al. 
(2004), normalised to the mean star-formation rate in CI 0024+16 and 
MS 0451-03 (again we note that these could be increased by ~ 18% to 
account for incompleteness in our color-selection of cluster members). 
The MS 045 1-03 point has an upper limit which corresponds to the 
extrapolated estimate for the total star-formation down to the luminos- 
ity limit in CI 0024+16 (6 x 10'° L©), while the data point shows the 
cluster's summed SFR down to the limit of our data (12 x lO'^L©). 
There is evidence for an increasing rate of activity in more distant 
clusters, as traced through their mid-infrared populations. However, 
there is also clear evidence for a wide variation in activity in even mas- 
sive clusters at a similar epoch. This suggests that the mid-infrared 
populations are sensitive tracers of environmental changes within the 
clusters. The errors on the total star-formation rates for CI 0024+16 
and MS 045 1-03 are derived from boot-strap resampling of the mid- 
infrared distribution, whereas for Perseus and A 1689 the errors are 
from counting statistics. Note that there are systematic uncertainties in 
all of the estimates depending on the specific choice of SED in the con- 
version from mid-infrared luminosity to star formation rate (the plotted 
uncertainties also do not take into account of the errors in the cluster 
mass estimates). Upper limits are equivalent star-formation rates for 
one detected at our luminosity Umit of 6x 10'° L©, (and extrapolated 
to account for the coverage difference in the survey) in those clusters 
where no sources were detected above this limit. 

At low redshifts, we use the nearby rich cluster Perseus, ob- 
served with IRAS (Meusinger et al. 2000). We sum over the 



known cluster members with Lfir > 6 X 10'" Lq, and con- 
vert their luminosities to SFRs, giving a lower limit to the to- 
tal star-formation rate in the cluster of > 22 M© yr"', over an 
area equivalent to 10 degrees^, which in terms of physical cov- 
erage is similar to the 2 Mpc radius used in CI 0024+16 and 
MS 0451-03. We use the mass estimate of Ettori, De Grandi 
& Molendi (2002) of 3.1 x 1O'''M0 and extrapolate to 2Mpc 
to determine the normaUsed star-formation rate per mass to be 
~3M0yr-VlO"M0. 

Fadda et al (2000) and Due et al. (2002) observed A 1689 
(z = 0.18) at 15 fj,m with ISO detecting 16 cluster members 
within 0.5 Mpc of the core. With the mid-infrared to total- 
infrared conversion used in this work, we find two galax- 
ies (detected with the LW3 15/xm filter) above our luminos- 
ity limit. In order to estimate the star-formation expected out 
to a radius of 2 Mpc we assume a radial profile for the mid- 
infrared population similar to that found for CI 0024+16 and 
extrapolate the star-formation rate from within a radius of 0.5 
to 2 Mpc. The resulting star-formation rate is ~28OM0 yr"' 
within 2 Mpc of the core. We use the mass from King et al. 
(2002), correcting to 2 Mpc, yielding a mass-normahzed value 
of ~3OM0yr-VlO"M0. A 2218 is another rich cluster at a 
similar redshift to A 1689, with z = 0.175, however the mid- 
infrared activity in this structure is much lower than in A 1689. 
The 15/im observations of Biviano et al. (2004) found nine 
members within a similar radius (< 0.4 Mpc), but only one of 
these corresponds to a star-forming galaxy, and the median in- 
frared luminosity is only 6x10** L0, or just 0. 1 M0 yr~' . A 370 
is at a similar redshift to CI 0024+16, at z = 0.37, but also ap- 
pears deficient in LIRGs. Only one cluster member was de- 
tected at 15/im. We plot these points as upper limits using 
our luminosity limit, again extrapolating to account for radial 
coverage out to 2 Mpc using the shape of the radial profile of 
mid-infrared sources in CI 0024+16. We use the cluster masses 
from Pratt et al. (2005) and Girardi & Mezzetti (2001) for 
A 2218 and A 370 respectively, extrapolating to find the mass 
within 2 Mpc. The mass-normalised upper-limits for A 2218 is 
30 M0 yr-'/lO''* M0 and for A 370 is 9 M0 yr-'/lO'^ M0. 

J 1888.16 (CI 0054-27) is at z = 0.56, and was observed by 
Due et al. (2004) using ISOCAM at 15/Lim. Six mid-infrared 
sources were detected which are confirmed cluster members, 
whilst a further two have redshifts suggesting they are mem- 
bers of infalling groups at a slightly higher redshift. Us- 
ing the calibration of Chary & Elbaz (2001), Due et al. find 
that all of their detected members are LIRGs, with inferred 
Lfir > 1.3 x lO" L0, and with individual SFRs in the range 
2O-I2OM0 yr"'. We re-calibrate the SFR using the method 
set out in §3.3, and find a conservative lower limit to the to- 
tal cluster star-formation rate to be > 70 Mq yr~' . The mass 
of the cluster is determined by Girardi & Mezzetti (2001), 
and we extrapolate to 2 Mpc, giving a total star formation rate 
>7M0yr-'/lO'4M0. 

We present the results in Figure 6, in which we plot the sum 
of the SFRs within ~2 Mpc of the clusters' cores, normalized 
using the estimated total mass of the cluster (either based on 
lensing or X-ray models). There is strong evolution in the 
star-formation rates in the clusters out to high redshifts, but 
it is important to note that several clusters seem to have very 
low-levels of activity, below our luminosity limits. For exam- 
ple CI 0024+ 16 appears significantly more active than A 370, 
which is at an almost identical redshift. This might be point to 
differing envirormiental influences between clusters being the 
dominant influence on the star-formation histories of in-falling 
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galaxies. Nevertheless, assuming that we have only selected 
those clusters with significant activity in Figure 6, there still ap- 
pears to be strong evolution in the sample out to z ~ 0.5. This 
observation could simply reflect the increase in number of star- 
forming galaxies seen in clusters and the field out to this red- 
shift - these obscured systems representing the high luminosity 
tail of the general 'blue' population. The idea that the clus- 
ter (obscured) star-formation histories mimic that of the field 
is supported by the rough consistency with the (1+z)^ trend 
found by Cowie et al. (2004) for the number of star-forming 
ULIRG radio sources out to z ~ 1.5. Kodama et al. (2004) 
present a similar analysis for the evolution of the Ha-derived 
star-formation in clusters out to z ~ 1. They find relatively 
strong evolution in the total-SFRs in clusters over this range, 
^ (1 +z)*, but as in this study, there is considerable scatter in 
the total star-formations rates between clusters, even after mass- 
normalisation. This hints that, although there might be a nom- 
inal rise in the level of star-formation in clusters out to z ~ 1, 
this is mitigated by the fact that individual cluster environments 
have a strong influence on the star-formation histories of their 
constituent galaxies. This scatter is due to the complexity of 
processes operating solely in the dense envirormients. A de- 
tailed study of the sources in CI 0024-H16 and MS 045 1-03 may 
elucidate this issue. 

5. CONCLUSIONS & SUMMARY 

We have used the MIPS instrument on Spitzer to survey 
the 24y^m populations of two optically rich clusters at z ~ 
0.5: CI 0024-1-16 and MS 0451-03. The samples are ~ 80% 
complete at 200/Jy, corresponding to total (8-1000/xm) in- 
frared luminosities of6x 10'" Lq and 12 x lO'^L© at z = 0.39 
and z = 0.55 respectively, equivalent to minimum SFRs of 
^ lOMgyr"' and ~ 2OM0yr-'. We detect a total of 986 
and 1018 mid-infrared sources above this flux limit across 
~ 25' X 25' fields in C10024h-16 and MS 0451-03. Our ob- 
servations probe from around the cores out to the turn-around 
radius at ~ 5 Mpc where the clusters merge into the field. In 
MS 045 1-03, we also analyse archival MIPS observations of 
the central ^ 5' x 5' of the cluster, which our observations had 
to avoid. Similarly in CI 0024-1-16 we make use of an ISO 15/xm 
survey from Coia et al. (2005) in the central region to build up 
a picture of the distribution of mid-infrared sources over the 
complete range of cluster environments. 

We exploit optical-near-infrared colors for the mid-infrared 
sources to reduce the background field contamination. We 
find a statistical excess of mid-infrared sources (within ~ 
5 Mpc of the cluster core) at 524^111 > 200/iJy associated with 
C10024H-16: 155±18. In contrast MS 0451-03 has a less sig- 
nificant population of mid-infrared sources, 28 ± 17, although 
we note that there are a small number of confirmed 24/xm mem- 
bers in MS 0451-03 in our on-going spectroscopic survey of 
this cluster. 

Using our deep optical and near-infrared imaging of both 
clusters we show that the 24/im sources in C10024H-16 are 
mostly associated with star-forming galaxies, with typically 
blue (B-R) colors, but which can be dust reddened by up to 
Ay ~ 2.5 mags. We also compare the infrared star-formation 



rate to that derived from an optical narrow-band Ha survey 
of this cluster from Kodama et al. (2004). Typically the Ha- 
derived rates underestimate the extinction-free infrared rate by 
> 5x, suggesting significant obscuration of the activity in this 
cluster We find that the level of obscuration for these individ- 
ual cluster galaxies is comparable to that found for LIRGs in the 
field at similar epochs. This suggests that starbursts in clusters 
are similar (at least in terms of extinction) as those triggered in 
low-density environments. However the variation in the 24/xm 
populations of C10024H-16 and MS 0451-03 suggests that the 
range of triggering and suppression mechanisms in clusters is 
complex. 

We estimate that the total star-formation rate (derived from 
the infrared) in the central region of CI 0024-1- 16 (< R200) is 
1000 ± 210 Moyr"'. MS 0451-03 is much poorer in mid- 
infrared sources, and we derive a total star-formation rate es- 
timate by summing over the smaU excess of objects, giving a 
total star-formation rate of 200 ± lOOM© yr~' within the same 
physical radius. We note however that our mid-infrared survey 
can miss some star-formation if a substantial number of lower- 
luminosity galaxies also exist in these clusters (as shown by the 
Ha survey of Kodama et al. 2004). However, we show that the 
majority of the activity is dominated by these dusty-starbursts. 

Finally, we look at the evolution of the specific star forma- 
tion rate per cluster with redshift from z ~ 0-0.5, using our 
new estimates for the total star-formation rates in C10024H-16 
and MS 0451-03. We compare these to estimates for lower and 
comparable redshift clusters studied with ISO and IRAS. We 
find that the high redshift clusters tend to have larger total star- 
formation rates compared to the more quiescent low-redshift 
ones, with an evolution similar to that of field LIRGs. However 
there is considerable scatter in this relation, and the evolution 
may only apply to the most active clusters. Although it is still 
unclear exactly what processes govern the star-formation his- 
tories of rich clusters, this study has shown that rich environ- 
ments can sustain significant amounts of hidden star-formation, 
and that this seems to increase at least out to z ~ 0.5. This hid- 
den activity may have a profound influence on the life-cycle 
of galaxies in high-density regions and the formation of the 
passive galaxy populations, ellipticals and SOs, which inhabit 
these environments at the present-day. We will investigate 
these issues in more detail in our next paper where we bring 
together spectroscopic and morphological information on the 
24/im population in these fields. We are also extending our sur- 
vey with new panoramic observations of distant clusters with 
Spitzer Cycle 3 GO time using both MIPS and the IRS. 
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Spitzer Science Center. J. E.G. is supported by a UK Particle 
Physics and Astronomy Research Council studentship; I.R.S. 
and G.P.S. acknowledge support from the Royal Society. 
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Spitzer MIPS surveys of two distant clusters 



Here we provide a brief description of the technique which we have adopted for identifying the optical counterparts to the mid- 
infrared sources in our catalog. We have chosen to apply the technique of de Ruitter, Willis & Arp (1977) who use a Bayesian 
estimator for the probabiUty, p(id | r), to identify whether an optical source in close proximity to a mid-infrared source is Ukely to be 
a true match, rather than a chance unrelated object. 

First we parameterize the radial distances of sources in terms of a dimensionless variable: 

r= + — (Al) 

yaa crs J 

where Aa{5) are the positional offsets of the sources in R.A. and Dec. (in the sense mid-infrared - optical), and (Taxs) are the 
uncertainties in the positions (we assume that the optical uncertainties are negligible). A likelihood ratio LR can be constructed 
representing the probability of finding a genuine match (id) compared to the probability of finding a confusing (c) source within a 
radius r. 



, P(r id) 1 
p(r I c) 2A 



r^2 

y(2A-l) 



(A2) 



where A = iraads'S, with S representing the surface density of optical sources down to some limit. The physical interpretation 
of equation A2 is that we are searching for a trade-off between finding a confusing background source, the distribution of which 
is governed by Poisson statistics; and the probability that the genuine match is located within dr of the reference source - this is 
described by the Rayleigh distribution. 

In our case A ^ 20E, and since we know that a cluster exists and the surface density of sources will not be uniform across the 
frame, we calculate the surface density of /?-band sources brighter than the nearest-neighbour - S(ot > wnn) - in annular bins of 
width 1' centered on the cluster at the radial distance of each mid-infrared source. The final part of the calculation requires us to 
calculate p{id \ r) - having found a source at some radius r, what is the probabiUty it is the genuine match? Bayes' theorem provides 
a way to estimate this value, but first we require an a priori estimate - namely what fraction 6 of mid-infrared sources have optical 
counterparts? To estimate 6 we count the number of sources with LR > 1, over the frame, denoting these 'matches' and compare 
to the number with LR < 1. 9 is insensitive to the choice of likelihood threshold, and is approximately 0.77 for C10024H-16 and 
MS 0451-03. We can then find: 

p(id|0==^^^ (A3) 

where X = 6/1-6. We define a positive match of sources at some cut-off value L of the likelihood ratio, which we optimize to 
provide the best completeness C and reliability R of the sample: 

LRi<L 



R = l-J2 Pi(c\r)/Nia (A5) 

LRi>L 

with A^id = J2Pi^^ I r)- We choose a value of p{id \ r) = 0.82 for positive matches, corresponding to /? = C ~ 97 % in both clusters. 
To improve the matching algorithm, we perform an identical calculation for the next-nearest neighbour to the mid-infrared sources. 
Thus for each MIPS source we have a matching probability for both its nearest and next-nearest optical neighbours. This allows 
us to flag possible mergers (where both probabilities exceed our match threshold). If pnnn > Pnn then we chose the next-nearest 
neighbour as the genuine match. 
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